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Abstract. - We study theoretically the capillary-gravity waves created at the water-air interface 
by an external surface pressure distribution symmetrical about a point and moving at constant 
velocity along a linear trajectory. Within the framework of linear wave theory and assuming the 
fluid to be inviscid, we calculate the wave resistance experienced by the perturbation as a function 
of its size (compared to the capillary length). In particular, we analyze how the amplitude of the 
jump occurring at the minimum phase speed Cmin = {Agj/ py^^ depends on the size of the pressure 
distribution (p is the liquid density, 7 is the water-air surface tension, and g is the acceleration due 
to gravity). We also show how for pressure distributions broader than a few capillary lengths, the 
result obtained by Havelock for the wave resistance in the particular case of pure gravity waves 
(i.e., 7 = 0) is progressively recovered. 



Introduction. — Water waves are both captivating 
and of great practical significance ^|-^- They have thus 
attracted the attention of scientists and engineers for many 
decades 4 . Water waves can, for instance, be generated 
by the wind blowing over the ocean, by a moving ship on 
a calm lake, or simply by throwing a stone into a pond. 
Their propagation at the surface of water is driven by a 
balance between the liquid inertia and its tendency, under 
the action of gravity or of surface tension (or a combi- 
nation of both in the case of capillary- gravity waves) ^ to 
return to a state of stable equilibrium 5 . The dispersive 
nature of water waves is responsible for the complicated 
wave pattern generated at the free surface of a still liquid 
by a moving disturbance such as a partially immersed ob- 
ject (e.g. a boat or an insect) or an external surface pres- 
sure source |6 . The propagating waves generated by the 
moving disturbance continuously remove energy to infin- 
ity. Consequently, the disturbance will experience a drag, 
called the wave resistance. In the case of ships (for 
which surface tension is negligible), this drag is known to 
be a major source of resistance 7 and has been analyzed 
in detail by Havelock \S^. The case of objects that are 
small compared to the capillary length has been considered 
only recently [9^16] and has attracted strong interest in 
the context of insect locomotion on water surfaces 17-19 . 
For such objects, one has to take into account both grav- 
ity and surface tension. In the case of a point-like surface 
pressure distribution, this leads to a discontinuity of the 



wave resistance at a critical velocity given by the mini- 
mum of the wave velocity Cmin = (4^p/7)^^^ 9 • For clean 
water at room temperature, one has Cmin ~ 0.23 ms~^. 
When the velocity V of the surface pressure distribution 
is smaller than Cmin, no steady waves are generated and 
the wave resistance vanishes. Emission of steady waves be- 
comes possible only when V > Cmin, leading to the onset 
of a finite wave drag 20 . The wave resistance disconti- 
nuity at the critical velocity Cmin has been experimentally 



investigated by several groups 13 , 14 



It is important to notice that while the analysis of ^ 
was mainly concerned with a point-like surface pressure 
distribution ^21j, real perturbations (like insects) have fi- 
nite sizes. The aim of the present paper is thus to analyze 
in detail the role played by the finite size of the surface 
pressure distribution for the wave resistance f22'. We will 
see in particular that for b much larger than a^"^, the in- 
fluence of surface tension is negligible. In that case, the 
results of Havelock [8^ for pure gravity waves are recov- 
ered. 

Formulation. — Consider an incompressible, inviscid, 
infinitely deep liquid whose free surface is unlimited. We 
take the x?/-plane as the equilibrium surface of the fluid 
and the z-axis along the upward direction perpendicular 
to the equilibrium surface. We study the wave motion 
created by an external surface pressure distribution that 
moves with speed V in the negative x-direction. In the 



p-i 



M. Benzaquen et al 



frame of that moving disturbance, physical quantities are 
stationary: the pressure distribution is given by P(x, y) 
and the displacement ( of the free surface from its equi- 
librium position is of the form ({x^y). 

In order to calculate the wave resistance experienced 
by the disturbance, we use a method first introduced by 
Havelock [s]. According to this author, we may imagine 
a rigid cover fitting the surface everywhere. The assigned 
pressure system P{x^y) is applied to the liquid by means 
of this cover; hence the wave resistance is simply the total 
resolved pressure in the x direction. This leads to 



R 



dx dyP{x,y)d:,C{x,y)- 



(1) 



For the sake of simplicity, let us restrict ourselves to the 
case of a pressure system symmetrical around the origin so 

1/2 

that P{x^y) = g{r) with r = (x^ + y^) )• The Fourier 
transform Pikx^ky) is then a function only of k and can 
be written as P{kx^ ky) = G(/c), where 



G{k) = / drrg{r) / dO 
Jo Jo 

= 27r / dr rg{r)Jo{kr). 
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Jo denotes the Bessel function of the first kind of zeroth 
order. It has been shown by Raphael and De Gennes in 
[9] that in such a case the wave resistance R reduces to 
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where for V > Cmin, X is defined by cos x = Cmin/V. Equa- 
tion Q is an important result as it predicts how the wave 
resistance varies as a function of the velocity for any pres- 
sure distribution. 

Finite size effects on the wave resistance. — 

Equation Q was studied in detail in 9 in the particu- 
lar case of a point-like pressure distribution g{r) = pd{r). 
We here consider the case of a pressure distribution of 
finite size, b. For instance, we can assume the pressure 
distribution to be Gaussian 



9{r) 



Equation ([3| then becomes 



G{k) = pexp 



(6) 



(7) 



By inserting the Fourier transform of the pressure field 
G{k) in equation Q, we obtain the wave resistance as a 
function of V/c^in (see Figs {!] and |2|. 




Fig. 1: (Color online) Plot of the wave resistance R in units 
of p'^ Hi{7Tj)~'^ , as a function of V/cmin- The red (dashed) curve 
corresponds to a point-like pressure source g{r) = pS{r), the 
black one (solid line) corresponds to a Gaussian pressure field 
del of size 6 = 0.04 /^"\ 
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Fig. 2: (Color online) Plot of the wave resistance R in units 
of p'^ k{7tj)~'^ , as a function of V/cmin- The red (dashed) curve 
corresponds to a point-like pressure field g{r) = p6{r), the 
black one (solid line) corresponds to a gaussian pressure field 
of size b = 0.4 a^"^ (upper graph) and b = 0.55 a^"^ (lower 
graph). 

The red (dashed) in curve in Fig. [l] corresponds to a 
point-like pressure distribution g{r) = pS{r)^ the black 
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one (solid line) corresponds to a gaussian pressure field 
^ of size b = 0.04 We clearly observe, for both 

curves, the discontinuity (or jump) at V = Cmin that we 
discussed earlier. In Figjl] the black curve (solid line) 
presents a maximum at Knax ^ {pb) which is the 
first consequence of the finite size effects. This separates 
the behavior of the wave resistance in two regimes: below 
^max, R increases with the disturbance velocity whereas 
for V > Fmax, R decreases with V. 

In Figjl] the disturbance typical size b is much smaller 
than the capillary length When b becomes of the 

same order of magnitude than tv~^ (see Fig|2|, the am- 
plitude of the jump, denoted by Ar in the following, de- 
creases as shown in Figjs) When the typical size b gets 
close enough to the capillary length the maximum 

value of the wave resistance is obtained for T^ax = Cmin- 
Such a situation is depicted in the lower graph in Figj2] 
Figures [4] and [5] give quantitative information on this sit- 
uation: we show that in the particular case of a Gaus- 
sian pressure field it occurs for 0.5 ^ ^ ^ 1.65 
Note that the above interval is only valid for the Gaussian 
pressure field ([6|. If, for instance, one uses a Lorentzian 
pressure source p{27r)~^b{b'^ o o 

above interval becomes 0.3 



ered. The amplitude of the jump Ar of the wave resistance 
in units of Hi{7rj)~^ , is depicted in Figjsjas a function of 
bhz. We notably observe that Ar is significantly suppressed 
when the typical size b becomes greater than 2.5 In 
such a case, the jump in the wave resistance might be 
difficult to observe experimentally (see also Fig Jg] below). 

One may also wonder how the maximum in the wave 
resistance scales with b. The abscissa (Knax) and ordinate 
(^max) of the maximum of wave resistance as a function 
of bhi are depicted in Fig j4] and FigjS] 



^)-^/2 instead of ([6|, the 
-1 <b<2.2K-\ 





Fig. 4: Plot of Vmax corresponding to the maximum of wave 
resistance in units of V/cmin, as a function of bn. The black 
dotted line (horizontal) corresponds to a situation in which the 
maximum value is reached at Vmax = Cmin- For ^ 1, Vmax 
scales as (pb), whether for bn^ 1, Vmax scales as 
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Fig. 3: Plot of the amplitude of the jump Ar of the wave 
resistance in units of p'^ Hi{7Tj)~'^ , as a function of bn. 
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From equation Q on can get the jump amplitude at 
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Equation ([8| is an important result as it gives how the 
jump of the wave resistance at the critical velocity Cmin 
varies with the size of the pressure distribution. In the par- 
ticular case of the Gaussian pressure field given in equation 
([g]), the jump amplitude becomes 



p- 



exp {-b^k^^ 



(9) 



,7r7 / 2^2 

Note that in the limit 6^0, the result of fo] A^i = 
7r/(2>/2) for of a point-like pressure distribution is recov 
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Fig. 5: Plot of the maximum of wave resistance i?max in units of 
p"^ Hi(7Tj)~'^ , as a function of bn. The dotted line corresponds to 
a situation in which the maximum value is reached at Knax = 



In both figures the dotted line corresponds to the situ- 
ation in which the maximum value is reached at l^ax = 
Cmin (see lower graph in Figj2|. 

In the limit bhc ^ 1^ the asymptotic behavior Fmax — 
^/gb can be simply recovered by substituting b~^ to k in 
the pure gravity wave dispersion relation uj{k)/k = ^J~gjk. 
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In the opposite limit <C 1, the asymptotic behavior 
Knax — yl/ipb) can analogously be obtained by sub- 
stituting to k in the pure capillary wave dispersion 
relation uj(k)/k = ^J^kj p [23^ . 

We observe in Figjs] that i?max (and hence the whole 
wave resistance) strongly decreases with the size h of the 
pressure source. Although this might be surprising at first 
sight, one has to notice that we have kept constant the 
magnitude p = J dx dy P(x, y) of the pressure source while 
varying b. 




Fig. 6: Plot of the quotient An/Rmax as a function of bK. The 
dotted line corresponds to a situation in which the maximum 
value is reached at Knax = Cmin- 

In order to discuss the general shape of the graph of the 
wave resistance, and whether or not the jump in the wave 
resistance can be easily detected, we plot the ratio between 
the amplitude Aji and the maximum of wave resistance 
^max (see Figj6|. One can see how the ratio A^^/i^max 
starts at zero (for 0, An is finite and i^max ^ oo) 

and increases until reaching the constant value 1 in the 
range 0.5 ^ ^ ^ 1.65 discussed earlier (see Figj2|. 
For larger values of 6, the ratio An/R^g,^ is essentially 
suppressed. 

Pure gravity waves. — As we have seen above, Vmax 
scales as ^/gb for bhz 1. One shall thus wonder if, more 
generally, our results for the wave resistance are for b 

well approximated by the result of Havelock |8| for 
pure gravity waves (7 = 0). 

If 7 = 0, the wave resistance Q reduces to 
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which, for the Gaussian pressure field (pj, yields 
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Fig. 7: (Color online) Plot of the wave resistance R in units 
of p'^ Hi(7Tj)~'^ , as a function of V/cmin- The black curve (solid 
line) corresponds to the formula for pure gravity waves (11), the 
red one (dashed) corresponds to the initial formula for capillary 
gravity waves (|4|. On the top we have 6 = 2 a^"^, in the middle 
we have b — 5/€~^, while at the bottom we have b — 10 k.~^ . 



Using equation (11), one can plot the wave resistance 
for pure gravity waves as a function of F/cmin for differ- 
ent values of b (see the black curves (solid line) in Figj?]). 
Note that Knax is now exactly given by \fgb. For veloci- 
ties smaller than 0.6\/^, the wave resistance is practically 
unnoticeable [24| . 

We have also plotted in FigjT] the wave resistance for 
capillary-gravity waves (red (dashed) curves). 

One can see how both curves come closer to each other 
as b is increased. For b much larger than ^ one does 
not see much differences between the two curves, meaning 
that in this limit the problem is essentially ruled by pure 
gravity waves theory. 



Concluding Remarks. — We have shown theoreti- 
cally how the finite size of an external axisymmetric pres- 
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sure source significantly modify the wave resistance, and 
in particular its singular behavior at the minimum phase 
speed Cmin, compared to the case of a point-like distur- 
bance [25] . Our study also provides a quantitative descrip- 
tion of the crossover between capillary-gravity and purely 
gravity wave resistance described in previous works. 



[23] Note however, 

[24] For velocities 
that equation 



that Vmax is not simply given by 



up to 0.8v^, 



on can show 
(11) is well approximated by 
l/x^)/(2x^) with X : 



//^(7r7-^)VW2exp(-l/x^)/(2x^) with x^V/ 
[25] In the context of animal locomotion on water surfaces, it 
would be interesting to generalize the results of the present 
study to non-axisymmetric pressure sources. 
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